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The hadronic quark structure is investigated in the frame of high energy electron 
proton scattering. A phenomenological model based on the Born approximation is 
used to calculate the transition matrix element for the quark system forming the 
proton target. A potential of electromagnetic nature is assumed for the calculation 
of the multiple scattering of electron with the constituent valance quarks of the pro- 
ton. It is found that the first two terms of the Born approximation are sufficient to 
describe the experimental data of differential cross section for the electron proton 
system only at low momentum transfer square. On the other hand, a two body scat- 
tering amplitude with a relevant form factor may give proper agreement extended to 
high momentum transfer square region. A harmonic oscillator wave function is used 
to reproduce the data of the experiment at low momentum transfer. This satisfies 
quietly the confinement condition of the quarks in the hadron system. However the 
data of the experiment at relatively high momentum transfer, show that hadronic 
quarks behave more freely and may be described by just a Coulomb wave function. 
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Experiments SLAC-E-140 = 1 - 3{GeV/cf = %GeV 

SLAC-E-136 g2 = 2.9 - 31.3(Gel^/c)2 Eg = 5 - 21.5Ge1/ 
SLAC-NH-11 g2 = 1.75 - 9>.9>'i{GeV/cf = 1.5 - %.9>GeV 

PACS numbers: 13.40.Fn , 25.40.Cm 
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I. INTRODUCTION 



The most recent data of the electron-proton (ep) scattering experiments at SLAC [1- 
3] have played a significant role in investigating the hadron structure. In the last decade, 
many trails have been executed to deal with the ep scattering problem. The Glauber multiple 
scattering approach [4] was used to present an eikonal picture for the ep scattering. The 
quantum chromodynamics QCD [5,6] was developed as the most convincing present theory 
for the interactions between those constituent quarks, formulated in terms of the exchange of 
colored vector gluons. At sufficiently high values of Q^, the running strong couphng constant 
is expected to become small enough, due to the properly of asymptotic freedom, to allow the 
use of perturbation theory to simplify QCD. However, there is considerable controversy as 
to how large a value of is sufficient for perturbative quantum chromodynamics PQCD [7] . 
The exclusive processes such as electron proton (e-p) scattering are predicted to have a single 
dimensional scaling [8] at large Q^. In this case only the valance quark states are important 
and a rough idea of the dependence can be gained by simply counting the number 
of quark-gluon vertices. Elastic form factor, for example, should scale asymptotically as 
(g2)-('^-i), where n is the number of valance quarks participating in the interaction. For (e- 
p) scattering n=3 thus the structure function behaves as q~^. Perturbative QCD [7] predicts 
calculable logarithmic departures from the dependence of exclusive amplitudes given 
by the simple dimensional scaling law. The earliest efforts [9] used unrealistic symmetric 
distribution amplitudes and required a large multiplicative factor to normalize the results to 
the data at ^ 10{GeV/cy. Chernyak and Zhitnitsky [10,11] proposed a model form for 
the nucleon distribution amplitude which satisfies the sum rules and in which the momentum 
balance of the valance quarks in the proton is quite asymmetric. The result of this model 
is justified in the diquark model [12,13]. The nonperturbative calculations succeeded in 
modeling the region < 20{GeV/cy fairly well but some difficulties have been found at 
large values of Q'^. By this article we aim to develop the eikonal optical picture of the (ep) 
scattering on the bases of the Born multiple scattering of the incident electron with the 
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constituent valance quarks of the target hadron. The paper is organized so that in section 
2, we present the hypothesis and the mathematical formalism of the model. Results and 
discussion as well the comparison with experimental data are given in section 3. 



II. THE MODEL 

During the electron-proton (e-p) scattering, we consider the proton as a bag including 
three valance quarks. We proceed using the semi-classical Born approximation which is a 
useful technique when the de Broglie wavelength of the incident particle is sufficiently short 
compared with the distance in which the potential varies appreciably. In this approximation 
we consider the multiple scattering series in which the projectile interacts repeatedly with 
the potential and propagates freely between two such interactions. We expect that the Born 
series converges if the incident particle is sufficiently fast so that it cannot interact many 
times with the potential and/or if the potential is weak enough. Starting with a plane wave 
as a zero order wave function of the incident electron, and by successive iteration, the nth 
order of the Born scattering amplitude is given by, 

/„ = -27r2<0fe,|C/|*„_i> (1) 

where fk, is the zero order plane wave and, 

*n = <f>kf{r) + J Go{r, f')C/(f)*„_i(f')rff' (2) 

and 

Go{f,f')^'-^^^^^^ (3) 

is the free Green's function. Let us denote by fg, ri, f2 and the coordinates of the incident 
electron and the constituent valance quarks of the target proton. According to the first Born 
approximation, the scattering amplitude is given by: 

/i = -(27r)^ / e'^-'''A{ro)dfo (4) 



/i = -(27r)2— / sin{qro)A{ro)dfo (5) 



This term represents the two-body scattering amphtude or the impulse approximation. The 
second Born approximation is, 

/2 = -(27r)2 / < kf\A\K > — — < n\A\ki > dn (6) 

which represents the double scattering term in the scattering amplitude. A(ro) is the average 
potential acting on the incident electron due to the target valance quarks ri with initial and 
final states and so that, 

Air,) = J ^b*{{ri})V{{r,})^a{{ri}){dfi} (7) 

In the present situation we consider that the electron-quark (eq) potential to be a pure 
Coulomb with the form, 

V{{n}) = ^ - ^ - ^ (8) 

roi ro2 Toz 

and Toi is the relative coordinate Voi = Vo — Ti The term l/voi may be expanded in terms 
of the spherical harmonics Yi^m as 

where r>(r<) is the greater (the lesser) of rg and rj. The radial integral of Eq.(3) may 
be divided into two parts: 

POO fTo roo 

- + (10) 

Jo Jo Jto 

(10) and substitute for r>(r<) by ri(ro) in the first integral and for r>(r<) by ro{ri) in the 
second integral. Then writing 

A{ro) = j:^A{n) (11) 
where A{ri) is the average potential due to the i^h quark, then 



A(ro) = / ^b*{{u})—^a{{u}){du} (12) 



Aiiro) 



f ^C{ri)—^a{ri)dfi 



= 47r/ \^{ri)\^-r^dri + 47r \^{ri)Wdri (13) 

^0 To Jro 

The evaluation of Ai{ro), {i = 1, 2and3} for the three quark system differs only in their 
signs due to the charge situation. On the other and the differential cross section is calculated 
either in terms of the multiple scattering Born approximation amplitude, 

da/dn=\f, + h + ...\^ (14) 

or in terms of the two body scattering amplitude with a proton form factor correction. 

da/dn = \fo\'F'{q) (15) 

here, fo is the scattering amphtude of the electron scattering by a point-like proton. The 
proton form factor F(q) is also an important physical quantity which reflects information 
about the particle structure. It is defined as: 

F{q) =< 0'|e^«"^"|C> > 

F{q) ^2t: j ^'{r)e''i'-^'^^^^{rydcos{e)dr (16) 

as the model is a phenomenological type, it is then branches in two ways, 

I) The first one assumes that quarks are confined inside the bag by an extremely deep 
potential, and the quark states are reasonably represented by a harmonic oscillator wave 
function. 

II) The second way assumes that the constituent quarks behave as almost free particles as 
is suggested by experiments of deep inelastic scattering [1,14]. However, in all cases the 
{e~p) scattering is treated as the collision between the high energy incident electron with 
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a composite proton system. Proceeding with the first assumption ( I ) and assuming a 
harmonic oscillator wave function for the valance quarks, 

*a = *6 = ( A)3/4,-.Va^ (17) 

Then the average scattering potential is; 

where r(3/2, 0, r^/a^) is the incomplete Gamma function of order 3/2. The scattering am- 
plitude and the differential cross section are to be calculated numerically. Moreover, the 
form factor is calculated for the proton system as, 

Fhosciq) = l/exp(a^qysf (19) 

Considering now the second case (II) of elastic scattering and that all quarks are in the IS 
ground state of a Coulomb wave function, 

= -^e-'/^ (20) 

'a' is the proton radius , then 



1_ ^ _a + 2r„ a + 2ro + 2r^ 



Consequently, 



fM - [f - - -atn{a g/2)] (22) 

da/dn = 4[^ - ^-p^ - iatn(a g/2)]^ (23) 
And the form factor is then, 

i^Co«K?) = 16/(4 + aV)' (24) 



7 



III. RESULTS AND DISCUSSION 



Data used in this article are those from the experiments coded SLAC-E-140 [1] at low- 
momentum transfer extended from 1 to 3 {GeV/c)'^ using electron beam of energy 8 GeV, 
and SLAC-E-136 [2] at a wide range extended from 2.9 to 31.3 [GeV/cY. The later was 
conducted using accelerated electrons with energies from 5 to 21.5 GeV were elasticaly 
scattered by protons in a liquid -hydrogen target at Stanford Linear Accelerator Center. 
The center of mass energy of the reaction is found in terms of the energies {Ef.^Ep) of the 
scattered electron and proton respectively and their opening angle 9. 



The differential cross section is calculated in the frame of multiple scattering using 
Eq.(14) for both the experiments SLAC-E-136 and SLAC-E-140. The result of calculation as 
well as the experimental data are displayed in Figs.(l) and (2) respectively. It is found that 
only the first two terms of the Born series give appreciable contribution to the differential 
cross section. The second term contributes for not more than 10 % in all cases. The cal- 
culations are carried out using harmonic oscillator wave function (solid lines) and Coulomb 
wave function (dashed lines) for the valance quarks forming the proton system. The exper- 
imental data are reproduced to a limited extend by the harmonic oscillator wave function 
only at low momentum transfer region. On the other hand the use of the Coulomb wave 
function shows better representation of the data at high momentum transfer. This result 
reflects the fact that quarks inside the hadron are quite confined at low momentum transfer 
region while at high momentum transfer, quarks behave more freely and the confinement 
condition becomes no longer necessary. However, the comparison with the experiment shows 
also that the multiple scattering approach couldn't give satisfactory agreement particularly 
at high momentum transfer regions. The differential cross section is recalculated in terms 
of the proton form factor according to Eqs. (15-16). Proceeding the same analogy, so that 
the proton form factor is calculated twice using the harmonic oscillator and the Coulomb 




(25) 
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wave functions. The result are displayed in Fig. (3) compared with data of the experiment 
SLAC-NH-11 [3] as well as the prediction of the simple dipole model of the following form, 

^ (l + gV0.71)2 (2^) 

The prediction of Eq. (15) is displayed with the experimental data in Figs. (4) and (5) for 
the experiments E-136 and E-140 respectively. It is clear now that the two body scattering 
amplitude corrected with a relevant form factor may reproduce the experimental data quite 
well, assuming harmonic oscillator wave function for the constituent quarks. In the above 
calculations, we use the units where {h = c = 1), hence the momentum in GeV/c may 
be represented in units of fm~^ with a conversion factor 0.2 In all the above cases the 
potential parameters are determined by the Chi- Square fitting method. A root mean 
square radius of the proton is found to be 0.8 and 1.0 fm on using the harmonic oscillator 
and the Coulomb wave functions respectively. The distribution of the hadronic system 
produced in e-p scattering is measured by ZEUS collaboration [15] at center of mass energy 
of 296 GeV. Comparison of the results with the QCD radiation has a strong influence on the 
characteristics of the final state. The data are reasonably reproduced by the Lund model 
based on a matrix element calculated in the first order of a, followed by a appropriate parton 
shower, as well as by the color dipole model. The HERWING parton shower model also gives 
a reasonable representation of the data. Neither the first order matrix elements alone nor 
the Lund parton shower model, without the matrix element calculation. 
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FIGURES 

FIG. 1. The differential cross section of ep scattering as calculated by the multiple scattering 

approach using for the proton valance quarks a harmonic oscillator wave function (solid line), and 
the Coulomb wave function (dashed line). The experimental data of SLAC-E-136 are represented 
by the cross signs (x) . 

FIG. 2. The differential cross section of ep scattering as calculated by the multiple scattering 
approach using for the proton valance quarks a harmonic oscillator wave function (solid line) , and 
the Coulomb wave function (dashed line). The experimental data of SLAC-E-140 are represented 
by the cross signs (x) . 

FIG. 3. The proton form factor calculated using for the proton valance quarks a harmonic 
oscillator wave function (solid line) , and the Coulomb wave function (dashed line) . The prediction 
of the simple dipole model is represented by (dotted line). The experimental data of SLAC-NE-11 
are represented by the cross signs (x). 

FIG. 4. The differential cross section of ep scattering calculated as two body scattering ampli- 
tude corrected with a relevant form factor using for the proton valance quarks a harmonic oscillator 
wave function (solid line), and the Coulomb wave function (dashed line). The experimental data 
of SLAC-E-136 are represented by the cross signs (x). 

FIG. 5. The differential cross section of ep scattering calculated as two body scattering ampli- 
tude corrected with a relevant form factor using for the proton valance quarks a harmonic oscillator 
wave function (solid line), and the Coulomb wave function (dashed line). The experimental data 
of SLAC-E-140 are represented by the cross signs (x). 
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